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Below the Curie temperature T c of a Heusler-alloy film, consisting of densely packed, but exchange- 
decoupled nanograins, the spontaneous magnetization M S (T) and static in-plane susceptibility x\\ (T) 
increase very slowly signalizing a suppression of magnetization fluctuations. The unpredicted vari- 
ation X|| (r) ~ G d (T), where Gd ~ Aff is the intergranular dipolar coupling, and also the magnetic 
freezing observed in the dynamic susceptibility at lower temperatures is quantitatively reproduced 
by Monte Carlo (MC) simulations with 10 4 dipolar-coupled moments on a disordered triangular 
lattice. At high temperatures, the MC spin configurations clearly reveal a dense gas of local vor- 
tex structures, which at low temperatures condense in regions with stronger disorder. This vortex 
depletion upon decreasing temperature seems to be responsible for the observed increase of the 
magnetic relaxation time. For weak disorder, the temperature dependence of the MC vorticity and 
a singularity of the specific heat at T v — \Gd{T v )/kB indicate a thermal transition from a vortex 
gas to a state with a single vortex center plus linear vortex structures. 

PACS numbers: 75.50.Tt; 75.40.Gb; 75.40.Mg; 75.70.Ak 



I. INTRODUCTION 

Thin ferromagnetic films are not only a very popular 
playground for physical research, but they also serve as 
the basic building blocks of the information technology 
(see e.g. Ref 1). The desire for high-density storage me- 
dia has motivated numerous studies on the preparation 
of nanomagnetic arrays and on the effects of the dipolar 
interaction between the nanomagnets on their dynamic 
behavior. In diluted arrays, the superparamagnetic fluc- 
tuations across the Neel^ energy barrier, E\, — KaV 9 , 
depending on the anisotropy Ka and the volume of the 
nanograin V g , limit the stability of the nanomagnetic mo- 
ments. Upon increasing density of the moments, the 
dipolar coupling between the grains, Ed, gains impor- 
tance and, due to its long range and angular anisotropy, 
the influence on the fluctuations and the dynamics of the 
nanomoments is extremely complex. This issue is of cur- 
rent interest and being studied from various directions. 

Recent theoretical work on the effect of structural 
disorder on the dipolar interaction in 2D ensembles 
of nanoparticles lead to the conclusion that the aver- 
age dipole energy decreases with increasing structural 
disorder 3,4 . On square/triangular periodic lattices^, the 
disorder destroys the antiferromagnetic/ferromagnetic 
dipolar ground state in favor of complicated non-collincar 
structures^ with complex energy landscape^ leading to 
the spin-glass like behavior. In a strongly coupled sys- 
tem with a large number of metastable minima an angu- 
lar rotation of a macroscopic magnetic moment of a given 
particle should cause neighboring particles to vary their 
magnetic directions as well, i.e. the cooperative rever- 
sal occurs^. While the static distribution of the energy 
barriers has been discussed in the literatureSiifliii, the 
description of the dynamic properties of a 2D interacting 
system is still lacking. Here, we mention the dependence 



of the freezing temperature on disorder, which is directly 
related to the temperature fluctuations. Moreover, crit- 
ical exponents in dipolar arrays and the dynamic sus- 
ceptibility have not been systematically studied. Con- 
sequently, it is yet not known whether in the presence 
of disorder the magnetic fluctuations are dominated by 
vortex-type excitations or via collinear reorientation of 
the magnetic moments. 

Upon increasing the coverage of a random assembly 
of nanometer-sized Fe-islands on CaFe2/ Si(lll) wafers, 
Scheinfein et al^2> observed a change from superpara- 
magnetism to in-plane ferromagnetism. Based on MC 
simulations, they attributed this long-range ordering 
to an interparticle exchange interaction, noting that in 
the absence of exchange such transition would be de- 
stroyed by thermal fluctuations. Much attention fo- 
cused on the hysteresis of the magnetization loops of 
two-dimensional (2D) arrays of nanomagnets, in par- 
ticular on the interplay between different kinds of par- 
ticle anisotropy, disorder, and the dipolar interaction. 
Hysteresis loops and coercive fields of 2D arrays of di- 
luted Co nanomagnets^ 3 ", of close-packed CoSiii, and of 
permalloy 8 particles were found to be in at least quali- 
tative agreement with numerical evaluations of suitably 
chosen models for structure and interactions. More re- 
cent numerical work considered smaller arrays of dipolar 
coupled particles with anisotropy on 2D lattices to obtain 
spin configurations, which display hysteresis loops with 
a multistep structure^. By ac susceptibility measure- 
ments in zero magnetic field, dipolar effects on the ther- 
mally activated magnetization dynamics have been de- 
tected on diluted hexagonal arrays of Co nanoclusters*^ 
and on close-packed, i.e. strongly interacting iron ox- 
ide nanoparticles 17 . In the first case, the slight enhance- 
ment of the relaxation time r and the energy barrier Ef, 
was found to be consistent with a model of independent 



anisotropic grains with weak dipolar interaction. In tl 
second case, stronger deviations from the Neel-behavi< 
have been interpreted in terms of a Vogel-Fulcher law fi 
r(T), i.e. spin-glass like collective dynamics^. 

Our work is devoted to the zero-field magnetizatk 
dynamics for a rather simple case that has not been pe 
sued in too great detail. We study a close-packed 2D a 
ray of magnetic moments, in which the dipolar couplii 
Gd between the moments is much larger than any oth 
energy, in particular, than volume, surface , and interfa< 
anisotropics. Experimentally, such system is realized 1 
a granular thin film, in which the grain boundaries c 
not transmit the ferromagnetic exchange. We invesl 
gate a thin film of the ferromagnetic cubic Heusler alk 
NiMnln, where the exchange between the magnetic M: 
moments arises from delocalized 3d-electrons provide 
by the almost nonmagnetic Nil&, which are also respo: 
sible for the potential of these alloys as spin-injectors 
lattice-matched semiconductors 1 ^ ^. Another nice fe 
ture of the Mn-based alloy is the small anisotropy fie 
(H a < 300 Oe)2£, most likely due to the weak couplii 
of the Mn-ion to the lattice. The large granular momer 
H g {T = 0) = 1.1 ■ 10 6 ^s, implies strong dipolar co 
pling, so that all experimental phenomena, in particul 
the magnetization dynamics are strongly collective. 1 
a unique feature of the present films, we note that, 1 
a variation of the temperature Gd(T) can be decreas< 
from the maximum value Gd(0)/ks = 3 • 10 4 K to ze: 
at Tq- A deeper understanding of the dipolar-dominatt 
magnetic fluctuations and of their dynamics is intend( 
by means of MC simulations^ of a closely related diso 
dered system. 

The outline of this work is the following. In Section 
the results of the magnetic measurements and of a stru 
tural characterization of the granular film are describe 
The main emphasis is devoted to the temperature varia- 
tion of the granular magnetic moment H g (T) and to the 
dc- and ac-susceptibilities in zero magnetic field, which 
probe the statics and dynamics of the magnetic fluctua- 
tions. By starting at Tq, the magnitude of the dipolar 
coupling, Gd{T) ~ fig(T), is increased up to large values 
to study its effect on the fluctuations. In Section III the 
MC simulations as a function of the normalized temper- 
ature T — ksT/Gd are presented. The results for the 
static and the time dependent susceptibility for different 
degrees of disorder are directly compared to the exper- 
imental data. Based on the surprising good agreement 
with respect to temperature variation and disorder, we 
then discuss the largely unexpected experimental find- 
ings in terms of the spin configurations provided by the 
MC simulations. Section IV contains the conclusions and 
some open questions. 



II. EXPERIMENTAL RESULTS 

The present film has been deposited on p-Si(lOO) by 
evaporating Ni and the stoichiometric Mnln alloy from 
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Figure 1: (a) TEM image of the NiMnln film with trian- 
gulation of the grain centers, (b) Distribution of grain areas 
obtained from the triangulation in (a) and of the grain bound- 
ary thicknesses (upper and right-hand scales, explained in the 
text). Solid curves are Gaussians, the dashed curve is a log- 
normal distribution with the same mean area A g = 680 nm 2 . 
(c) Temperature dependence of the resistance of the NiMnln 
film on p-Si (from Ref 26). 



two independent sources as described in more detail 
elsewhere^. X-ray spectroscopy and diffraction revealed 
the so-called half-Heusler alloy NiMnln, which can be 
described as the L2\ structure of the full-Heusler alloy 
Ni2MnIn (see e.g. Ref. l25lh where one of the four inter- 
vening fee lattices for each metal ion is missing. The 
granular structure of the film has been evidenced by 
atomic force microscopy (AFM) and also by transmiss- 
sion electron microscopy (TEM) on a film simultaneously 
prepared on amorphous carbon, see Fig. 1(a). AFM scans 
displayed very sharp film edges and provided for the 
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Figure 2: Temperature dependence of the technically sat- 
urated magnetization M a of a NiMnln Heusler-alloy film. 
M S (T) is defined by the plateaus of the magnetization 
isotherms below the Curie-temperature, shown by the in- 
set. The solid curve fits the data to the power law, M a ~ 
(f - TfTofi with the lar S e exponent 0=0.66. 



thickness t=35±l nm. We analyzed the TEM image by 
means of a triangulation : taking the grain centers , i.e. 
the vertices shown in Fig. 1(a), as lattice points, one finds 
the distribution of the intergranular distances d, which 
define the areas of the hexagonal grains, A = d 2 \fi/2. 
The results depicted in Fig. 1(b) are consistent with both 
the gaussian and the log-normal distributions and the 
same mean grain area A g = 680 ± 40 nm 2 . This im- 
plies for the mean distance between the grain centers 
d g = 28.0 nm. We will refer to these results below, when 
analyzing th ac-susceptibilities. 

Using a digital processing of the TEM image, the con- 
trasts produced by the grain boundaries have been deter- 
mined and analyzed by a Fast Fourier Transformation 2 ^. 
The results shown in Fig. 1(b), display a rather sharp 
Gaussian peak at the spatial frequency 2n/A = 
5.77 nm~ 1 , which implies for the mean grain boundary 
thickness A = (l.l±0.1)nm. We believe that this rather 
large value for A is responsible for the negligible contri- 
bution of the Heusler film to the resistance curve R(T) 
presented in Fig. 1(c): the plateau of R at low tempera- 
tures and the drop near room temperature arise fom a 
highly doped p + — Si layer and electron-hole excitations 
in the Si-wafer, respectiveljiSi. The accuracy of the data 
sets a lower (nonmetallic) limit of p > 1 mflcm to the 
Heusler-film. As the ferromagnetic exchange between lo- 
calized M n-moments in the cubic structure results from 
the conduction electrons^, such coupling across the wide 
and disordered boundaries appears to be rather unlikely. 

The magnetic measurements on the NiMnln film 



(11 mm 2 area) have been performed by a SQUID magne- 
tometer (MPMS2, Quantum Design, San Diego). Imme- 
diately below the Curie temperature of NiMnln, the mag- 
netization isotherms rise rather steeply to their technical 
saturation value M S (T) (see inset to Fig. 2). According 
to Fig. 2, the temperature variation of M s obeys rather 
accurately the power law M S (T) = M a (0)(l - T/T c f 
with M s (0) = 5.1 kOe. As the most striking result we 
note that the exponent (3 = 0.66 ± 0.02 is not only sig- 
nificantly greater than (3 — 0.35 for the 3D Heisenberg 
ferromagneljSi but also larger than the mean field value 
j3 = 0.5. We conjecture here an effect of the fluctuating 
dipolar fields in the disordered granular structure, but 
cannot yet offer a more detailed discussion of this inter- 
esting observation. 

Basically, fluctuation effects should be more pro- 
nounced in the initial susceptibility, the temperature de- 
pendence of which was determined more accurately by 
ac-measurements and the zero-field cooled (ZFC) mag- 
netization, see Fig. 3 (a). Upon decreasing temperature, 
the polar ('out-of-plane') susceptibility reaches soon the 
demagnetizing limit x±. — lj whereas the parallel ('in- 
plane') one, X\\t increases up to frequency dependent 
peak temperatures Tf. The enveloppe on the high- 
temperature side of the peak defines the static limit 
Xy(T) = x\\{T,u> — > 0), which is analyzed by a Kouvel- 
Fishcr plot™ in the inset. The appearance of a straight 
line for [d(lnx||)/dT] signals a power law, X\\(T) = 
X||(0)(1 - T/TcY, the Curie temperature T c = 323.5 K 
and exponent k = 2.6±0.1 of which are given by the 
intercept with the T-axis and by the slope, respectively. 
Similar as for M S (T), the large exponent k characterizes 
a rather slow increase of the susceptibility, as compared 
e.g. to the superparamagnetic limit, XII p 1 ~^ ^c) ~ 
/i 2 (T) - (1 - T/T c ) 213 , with 2f3 = 0.70 expected for non- 
interacting granular moments ^ g (T) — M s (T)V g , where 
V g is the mean grain volume. This feature indicates a 
strong suppression of magnetic dipole fluctuations by the 
intergranular dipolar coupling. 

In Fig. 3(b), the freezing temperatures T'j and T'j fol- 
lowing from the peaks of the real (Fig. 3(a)) and imag- 
inary (Fig. 4) parts of X\\(T, u = 27r/), respectively, are 
analyzed in terms of the traditional Neel-Arrhenius plot 
for a thermally activated relaxation time of the magne- 
tization, t(T) = To ex-p(Ef,(T)/kBT). Here the usual 
1 /T variable has been modified by taking into account 
a temperature dependence of the energy barrier which 
is demanded by the fact that otherwise a completely 
unphysical attempt frequency t^ 1 = 10 146 s _1 is ob- 
tained. We argue that the barrier is related to the dom- 
inant anisotropy in the granular film, i.e. to the dipole- 
dipole interaction between the granular moments, be- 
cause, as mentioned previously, the crystalline anisotropy 
field of the Heusler alloy is much smaller than the spon- 
taneous magnetization M s . Hence, we assume Eb(T) ~ 
/i 2 (T), for which the temperature variation of M S (T) = 
Pg(T)/V g in Fig.2 implies E b (T) = E b (0)(l - T/T c ) 20 . 
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Figure 3: (a) Temperature variation of the real part of the ac- 
susceptibility x'\\ or NiMnln measured in zero magnetic field 
between 0.03 Hz and 928 Hz for ac-fields parallel to the film- 
plane. Also included are the zero-field-cooled in-plane suscep- 
tibility, Mzfc/H, recorded during warming in H=l Oe and 
the out-of-plane susceptibility f° r 10 Hz. All solid curves 
are drawn as guides to the eye. Inset: Kouvel-Fisher analysis 
of the static in-plane susceptibility, X\\(T) = x'\\{T,uj — > 0). 
The straight line defines the Curie-temperature Tc and the 
exponent k=2.6(1) of the power law x\\{ T ) ~ (1 - T/T C ) K , 
depicted by the broken curve in the main frame, (b) Arrhe- 
nius analyses of the peak temperatures Tf of x'\\ ( d marked 
by arrows in panel a) and x" C m Fig- 4) for a temperature 
dependent energy barrier, E b (T) = E b (0)(l - T/T c ) 2/5 . 



Within the Debye model, 



X\\(u,T) = X\\(T)— —. jTpr, 

' " 1 + IUJT(1 ) 



(1) 
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Figure 4: Temperature variation of the ac-absorption of the 
NiMnln film. The solid curves have been calculated from 
Eq.2, using a gaussian distribution of energy barriers of width 
a"6=0.24 and the mean barrier Eb(T) following from Fig. 3(b). 
The dashed curve has been calculated using the log-normal 
distribution depicted in Fig. 1(b). 



dependent MC simulations of the problem, to be pre- 
sented in section III, provide only the real part of the 
dynamic susceptibility. 

The present discussion of the dynamics in terms of a 
thermally activated process is also supported by the ocur- 
rence of a blocking of the ZFC magnetization at 7], in 
Fig. 3(a). Depicting T b on the Arrhenius plot for Tj. 
in Fig. 3(b), a blocking time T b = (2Tif b )~ 1 = 24 s is 
obtained, a magnitude which is usually accepted for the 
onset of magnetic irreversibility^.. A detectable coercive 
field, H c = 2 Oe, appears near 270 K, i.e. below the 
temperature range of interest in this work. 

The temperature dependence of the absorption compo- 
nent of x\\ is exemplified for three frequencies in Fig. 4. 
As the width of the peak turns out to be much larger than 
for the pure Debye-shape, Eq.l, a distribution of energy 
barriers in the disordered granular structure is very likely. 
Based on the results in Fig. 1(b), we assume a gaus- 
sian distribution of granular blocking volumes V about 
a mean value V g , P(v) — (2na v )~^ exp(— (v — l) 2 /2a v ), 
where v = V/V g . This implies for the barrier Eb(v) = 
v- Ei,, where E b is the mean barrier of the grains. The ac- 
susceptibility computed with t(T, v) — tq exp(vE b /ksT) 
and a v = 0.24 from the analysis in Fig. 1(b), 



the peak temperatures of x' an d x" obey the relation 
E b (Tf)/T f = ln(/o//). Hence, the slopes of both plots 
in Fig. 3(b) determine the barrier at zero temperature, 
Eb(0)/k B — 150 • lO 3 -^, while the displacement of both 
lines arises from the difference /o//o ~ |l nwr o| ~ 3.5. 
The preexponential time tq = 1.4 • 10~ 13 s follows di- 
rectly from the extrapolation of the x"-pea,ks to 1/ /q — 
2ttt exp(l/| lnwrol) = 27tt at 1/1/ = 0. We have out- 
lined this technical point in some detail, because our time 



X\\(u,T;a v ) = / dvP(v)v x\\(u,T,v) (2) 



agrees excellently with the data in Fig. 4. Hence, the 
results are fully consistent with the gaussian distribution. 
Alternatively, inserting the log-normal distribution from 
Fig. 1(b) to Eq.(2), we find also good agreement with 
the data, indicated by the dashed curve in Fig. 4, exept 
in the wing on the low temperatur side of the x"-peak. 
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These findings permit the estimates of two important 
quantities, i.e. the mean magnetostatic energy (in SI 
units fio — 4-7T • 1CP 7 Vs/Am) of a single grain, E g — 
^Ho M^Vg with the mean volume of a grain V g = A g t = 
(23.8 ± 1.5) 10 3 nm 3 , and of the mean dipolar coupling 
energy between adjacent grains, 



100 



G d (T) 



4-7T 



(3) 



At first, we compare the magnetostatic energy at T=0, 
E g /k B = 192 • 10 3 K with the energy barrier at T=0 de- 
termined from the Arrhenius analysis, E^{T — 0)/kg = 
150 • 10 3 K, to obtain E b = 0.78 E g . This result suggests 
a close relation between the barrier and the magneto- 
static energy, which for an out-of-plane flip of a granular 
moment sets an upper limit for E^. Of course, this up- 
per limit may not be reached due to a reduction of the 
in-plane dipolar field by the disorder and/or by thermal 
fluctuations of the adjacent moments. 



III. MONTE CARLO SIMULATIONS 

In order to explore the intriguing question to what ex- 
tent the disorder and thermal fluctuations affect the mag- 
netization statics and dynamics, we performed extensive 
MC simulations for 100 x 100 classical spins with dipolar 
interaction on a triangular lattice. The lattice sites have 
been randomized so that gaussian distributions of gran- 
ular distances d g with various widths aa were realized 
and, hence, the central features of the granular structure 
of the film were taken into account^ 

The interaction Hamiltonian reads 



where n 




, (Si • Tij) (Sj ■ Tij) 



= J d g is the mean distance between the mo- 
r, and Tj normalized to the nearest neighbor 



ments at 

spacing d g , and Gd denotes the strength of the coupling, 
given in Eq. 3. Our aim is to give a reasonable theoretical 
description of finite arrays. For that reason and in order 
to avoid symmetry adapted structures, open boundary 
conditions have been used. To prevent artificial effects we 
used no cutoff. A standard MC technique was used^ 3 .. In 
contrast to MC schemes for exchange coupled systems, 
where only the restricted movements of magnetic mo- 
ment are often used 3 ^, the rotational space was sampled 
uniformly , i.e. a moment can try any new angle. This 
is especially important in dipolar systems as the dipolar 
interaction often favors large angles between neighboring 
spins. An extremely slow annealing procedure with up 
to 150 temperature steps has been applied. The number 
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Figure 5: (a) In- and out-of-plane susceptibilities \\\ an d X-L 
evaluated from Monte Carlo (MC) simulations for dipolar cou- 
pled moments versus temperature normalized to the coupling 
constant, T — ksT '/Gd- The MC results are obtained after 
tMC = 5 ■ 10 3 MC steps for 100x100 classical spins on site 
disordered triangular lattices with gaussian widths ad = 0.15 
and 0.35 (thin curve). Above the peak temperature T p , the 
log-log plot reveals the power law X\\' C (T > %) = 4T~ 2 
(broken line), which is in absolute agreement with the exper- 
imental data for \\\ a t 30 mHz inferred from Fig. 3(a) using 
G d (T) = Gd{0)(l-T/T c y- 32 . (b) Arrhenius plot of the peak 
temperatures T p vs. MC time tuc yielding the energy barrier 
E¥ c for two different widths Gd- 



of MC steps (MCS) at each temperature, i.e. the MC ob- 
servation time tMC, has been varied between 0.5 • 10 3 and 
10 4 , so that the whole procedure contained up to 1.5 • 10 6 
MCS. To avoid metastable states we have performed two 
different simulations of the same system simultaneously, 
starting them at different seeds for the random number 
generator to ensure that the samples take different paths 
to the equilibrium. Only when both samples reached the 
same stable energy level it has been accepted that the 
system has reached equilibrium. 

The central results for the in-plane and perpendicular 
susceptibilities, X a° = 4nG d {0)/k B T ■ ((S 2 a ) - (S a f), 
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a = | and _L, are depicted in Fig. 5(a) versus the nor- 
malized temperature T — k B T '/G4. Very similar to the 
experimental observation in Fig. 3(a), X|| displays a 

peak at some temperature T p , the position and width 
of which depend on the strength of the disorder. On the 
high temperature side of the peak, the double logarith- 
mic plot clearly reveals a power law, xf ' (T > T p ) — 

4 T~ 2 . A variation of the MC observation time tuc 
does not change this behavior, so that this power law 
represents to the equilibrium susceptibility and, hence, 
a direct comparison with the experimental results for 
X\\ (T) can be made. For this purpose, we employ our 
x'-data for oj/2tt—3Q mHz from Fig. 3(a) and trans- 
form the temperature scale to the reduced MC-scale 
by means of Eq.(3), f = k B T/G d {0){l - T/T c ) 132 , 
where G d (0) = ^ M 2 (0)V 2 /And 3 g = 33.1 ■ 10 3 k B K. Us- 
ing this relation between T and T, we see immediately 
that the MC-result xf C (T > T p ) = AT' 2 implies di- 
rectly x|| (T) = x(0)(l - T/T c ) 2 - 64 . As a matter of fact, 
xf {T > T p ) agrees not only with the T-dependence of 
the static suceptibility in Fig. 3(a), but we find also per- 
fect agreement for the amplitude X|| (0) and for xx(T). 

Before discussing the origin behind this suppression of 
the magnetic fluctuations in the susceptibility by help of 
the MC spin configurations, we examine the maximum of 
xf i i- e - freezing behavior occuring at T p . To this end, 
we calculate xf C (T) for MC times between tuc = 5- 10 2 
and 10 4 MCS and find the peak temperature to decrease 
for longer tuc- The slopes of the Arrhenius plots for T~ x 
presented in Fig. 5(b) define the normalized energy bar- 
riers 

d{hitMc)/ dT p X = Eb IC /G d , which turn out to be 
reduced by increasing disorder ad- In order to compare 
these results with the experimental barrier, we note that 
ad measures the variance of the intergranular distances, 
so that the corresponding variance of areas is given by 
aA = 2&d- Hence, the MC barrier for ad = 0.15 may 
be compared to E b = (192 ± 10) 10 3 £;h K of the Heusler 
film with <7b — 0.24 extracted from the x"(T, w)-shape 
in Fig. 4, with which the width of the distribution func- 
tion P(v) in Fig. 1(b) was consistent. In fact the MC 
barrier for this disorder, E^ IC {0)/k B = 5.8 • G d {0)/k B = 
191 • 10 3 K, is almost identical with the measured value. 
However, regarding the slight difference between the dis- 
order strengths this identity appears to be accidental. 

Motivated by the overall agreement between the ex- 
perimental and the MC susceptibilities in the equilibrium 
and blocking regimes and for a deeper understanding of 
the observations it is rather suggestive to explore also 
the MC-generated spin configurations. Since the dipole 
interaction between the moments prefers the formation 
of vortices^, we introduce the socalled vorticity of the 
triangular lattice cell at r.;, 

Q z (r i ) = ^--rot z S(r i ). (5) 



Since at the temperatures of interest here the S z - 
fluctuations are negligible, the other components can be 
ignored, Q x = Q y = 0. 

At first, let us discuss the global, i.e. average of the 
vorticity modulus of the disc, 

N 2 

Q* = N- 2 J2\Qz(n)\- (6) 

i=l 

For vanishing disorder, the groundstate consists of a sin- 
gle vortex or antivortex localized at the center, Q z (v = 
0) = ±1 ( see Fig. 6(c)), hence Q z — A~ 2 turns out 
to be very small. At finite ad, Ed favors non-collincar 
spin configurations in regions with enhanced disorder. 
At low temperatures this is evidenced in Fig. 6(a) by 
the strong increase of Q z with rising ad- For the weakest 
disorder, the average vorticity Q z displays at T v — 0.5 a 
rather discontinuous transition upon increasing temper- 
ature. Within the accuracy of the data, this transition 
to a state with a maximum vorticity, Q z = 1, does not 
depend on the disorder. 

More specifically, the presence of a thermal transition 
is born out by the normalized specific heat Gy = ((Ej) — 
(E d ) 2 )/(k B T) 2 , see inset to Fig. 6(a) for a d = 0.05. 
Approaching T v from above, a power law singularity is 
found, Cy(T) - (f - f v )- a , with a = 1.73 ± 0.02, 
which is independent of tuc- This exponent is very 
large as compared to values for isotropic 3D ferromag- 
nets, a < 0.1—. Taking also the slow increase of the 
thermal susceptibility into account, it is suggestive to at- 
tribute this behavior to the vortex fluctuations at high 
temperatures. As the disorder increases, the peak of 
Gv(T) is observed to remain near T v = 0.5, but be- 
ing strongly suppressed. This indicates that the disor- 
der reduces the thermal fluctuations by freezing vortices 
at sites which are favored by the interplay between the 
structural disorder and the dipolar interaction. 

These effects of temperature and disorder on the vor- 
ticity Q z and the specific heat can be further illustrated 
by vorticity maps depicted in Fig. 6(b) for an interme- 
diate disorder, ad = 0.15, i.e. close to the experimental 
value. At low temperatures, a network of local vortex 
fluctuations is realized, the structure of which reflects 
the variation of the disorder assumed for the MC simu- 
lations. A more detailed inspection of the spin structure 
and the underlying lattice disorder reveals, that the nu- 
cleation of vortices is preferred at sites where the lattice 
spacing is reduced and that the dipolar interaction fa- 
vors a lining-up of vortices of the same sign. Such linear 
structures are clearly seen in Fig. 6(b) at the lower tem- 
peratures. As they separate vortex-free domains with 
almost homogeneous magnetizations pointing in reverse 
directions, these 'vortex-lines' act like domain walls in 
ferromagnets, and it is plausible to attribute the obser- 
vation of shorter relaxation times, i.e. smaller energy 
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Figure 6: (a) Temperature variations of the vorticity modu- 
lus, Eq. 6, indicating a thermal transition at T v — 0.5 for 
vanishing disorder ad- This feature is supported by the diver- 
gence of the specific heat shown by the inset for ad = 0.05. 
Vorticity maps determined from MC spin configurations at 
temperatures in the fluctuating (T ~ T p ) and in the frozen 
state (T <C T p ) for gaussian disorder strengths, (b) a"d=0.15, 
close to the experimental value, and (c) ad = 0. The diame- 
ter of the circles increases with the modulus of Q z (fi). The 
bright regions in panel (c) are free of vortices and display a 
spin circulation S indicated for T = 0.15. 



barriers, at larger disorder ad (see Fig. 5(b)) to the pres- 
ence of more walls. At the highest temperature, the en- 
tire disc is continuously covered by vortex fluctuations, 
which reduce further the relaxation time so that the equi- 
librium susceptibility, \\\ {T) = 4 T~ 2 , is reached in both, 
the experiment (Fig. 3(a)) and the MC simulation (Fig. 
5(a)). On the other hand, since these dipolar-induced 
vortex fluctuations are non-collinear spin structures, the 
susceptibility is strongly suppressed as compared to the 
superparamagnetic limit. 

Additional physical insight is provided by Fig. 6(c), 
where in the absence of disorder the thermal excitation 
of vortices is realized. Far below the transition at T v , i.e. 
in the groundstate, there is a single vortex in the center 
of the disc determining the direction of the circulation 
of the almost homogeneous magnetization M ~ S in the 
remainder of the sample. Just below T v ,atT — 0.45, the 
number of vortices is much smaller than in the presence of 



disorder (Fig. 6(b)). At the larger temperature, T = 0.9, 
the number of vortices increases dramatically resembling 
a gas of thermally excited vortices. We believe that an 
explanation of the vortex transition at T v = 0.5 should 
take into account the short range order of lining-up the 
vortices which appears at T = 0.45 in Fig. 6(c). On 
this basis, a Berezinskii-Kosterlitz-Thouless transition 35 
as conjectured by Prakash and Henley^i for 2D dipolar- 
coupled spins seems to be unlikely. We do not see any 
evidence for a binding of vortices and antivortices, but 
rather the trend of an attraction between vortices of one 
kind to form linear arrays. 



IV. SUMMARY AND CONCLUSIONS 

Upon decreasing the temperature below Tq of a thin 
Heusler alloy film consisting of exchange-decoupled fer- 
romagnetic nanograms, V g ~ (29 nm) 3 , the effect of 
the strongly increasing dipolar interaction between the 
grains on the ordering of the film has been investi- 
gated. The rather slow rise of the zero-field susceptibility, 
X\\(T) ~ Hg{T) ~ G 2 d (T), is in quantitative agreement 
with MC simulations of dipolar-coupled classical spins on 
a triangular lattice. We should emphasize, that no finite 
size effect was observed when reducing the system from 
100 x 100 to 60 x 60 spins. Regarding the comparison 
to the experiment, an effect of the circular sample shape, 
chosen for the simulations, is unlikely, since Fig. 6 re- 
veals, that at the temperatures of interest (T > T v ) ther- 
mal and disorder fluctuations of short range dominate the 
spin configurations. The MC spin configurations allow to 
attribute this rather striking feature to the thermal ex- 
citations of dense, small-scale vortices in the film plane, 
which are demanded by the dipolar coupling and which 
suppress the dipole fluctuations probed by the suscep- 
tibility. The temperature variation of the vorticity and 
of the specific heat, determined by the MC simulations 
in the limit of weak disorder, suggest the presence of a 
thermally induced phase transition at T v = \Gd{T v )/kB 1 
which appears to be driven by vortex fluctuations of yet 
unknown structure. 

The dynamic susceptibility of the granular film reveals 
a thermally activated magnetic relaxation across an en- 
ergy barrier E^. This barrier turns out to be slightly 
smaller than the mean magnetostatic energy of a grain, 
Eb < 2-xGd, and the time-dependent MC simulation re- 
produces this measured value Ej, also almost quantita- 
tively. This agreement has been achieved for a gaus- 
sian site disorder on the triangular lattice with a vari- 
ance close to that of the structural disorder of the film. 
This disorder of the film has been inferred from both the 
width of the measured ac-absorption and of the trian- 
gulation of the granular structure. As a matter of fact, 
the MC work predicts the energy barrier to be lowered 
by increasing disorder, which appears to be associated 
with the larger number of vortices generated by the dis- 
order. Due to the dipolar interaction, vortices tend to 
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form linear arrays, which act similar as domain walls in 
ferromagnets. Qualitatively, our observation is in anal- 
ogy to the increase of the barrier reported for diluted 2D 
arrays with increasing GdF^- 

The present results on the effects of thermally and 
structure induced vortex fluctuations on the statics and 
dynamics of the magnetization as well as the nature of 
the transition near T v call for further efforts from both 
sides, experiment and theory. Experimentally, a direct 
proof of the vortices could perhaps be provided by near- 
field magnetic force microscopy. Another challenge could 
be the dynamics of the dipolar-coupled disordered grains 
in the presence of a magnetic field, in particular, their 
hysteretic and switching behavior at low temperatures, 



T <C Tc- Finally, we should mention the yet unexplained 
reduction of the increase of the spontaneous magnetiza- 
tion, m s (T) ~ (Tc — T)' 3 , due to the large exponent 
(3 = 0.66, where an effect of the disorder is not unlikely 
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